Hamiltonian anomalies of bound states in QED. 
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The Bound State in QED is described in systematic way by means of nonlocal irre- 
ducible representations of the nonhomogeneous Poincare group and Dirac's method 
of quantization. As an example of application of this method we calculate trian- 



gle diagram Para — Positronium — > 77. We show that the Hamiltonian approach 



to Bound State in QED leads to anomaly-type contribution to creation of pair of 
parapositronium by two photon. 
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1. INTRODUCTION 



The bound states in gauge theories are usually considered in the framework of represen- 
tations of the homogeneous Lorentz group in one of the Lorentz- invariant gauges [1]. In this 



paper, we suggest a systematic scheme of the bound state generalization of S-matrix ele- 
ments, which is based on irreducible representations of the nonhomogeneous Poincare group 
in concordance with the first Quantum Electrodynamic (QED) quantization [2, 3] and first 
QED description of bound states [4, 5]. We obtain bound states by means of excluding time- 
component of four-potential [2, 3, 6] and Hubbard-Stratonovich transformation [7, 8], that 
is in agreement with general principles [9] . The time component is chosen in correspondence 
with the Markov- Yukava constraint of irreducibility [10, 11]. 

The aim of the article is to research the experimental consequences of such the Poincare 
group irreducible representations of QED (see also review [12]). 

As a test of this scheme we calculate process P — > 77, where P is parapositronium, that 
describes triangle diagram, with one real and the other virtual photons. 

The structure of the article is as follows. In section 2 we exclude time-component of 
four-potential from QED action. For derived action in section 3 we make the Hubbard- 
Stratonovich transformation and positronium field appear. Than we make semi-classical 
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quantization of the resulting system in section 4. In section 5 we calculate the triangle 
anomaly. And in section 6 we calculate the contribution in process 77 — > PP inspired by 
triangle anomaly. 

2. EXCLUDING TIME COMPONENT OF IN QED 

We start with usual Quantum Electrodynamic Lagrangian: 

Cqed = -\f^ - A^f - ^(iyfy - m)1>. (1) 

Below we will work in frame of reference in which the bound state as whole is at rest. 
The Lagrangian (1) contains nonphysical degree of freedom. This degree can be excluded 
by substitution the solution of classical equation of motion in action. 

Variation action providing by (1) over A leads to one component of Maxwell equation: 

AA - d d k A k = -j , 

where A = d±di + d 2 d 2 + <9 3 <9 3 . If we take the gauge: 

d k A k = (2) 

than the previous equation simplifies: 

AA = - Jo . 



The solution of this equation is: 



Ao = ~^Jo (3) 



where 

( x> y,z ^ = J dx ' d v' dz ' ^ x ' y ' z 



J(x - x') 2 + {y- y') 2 + {z- z') 2 

K 3 \{(x,y,z)} 

Under gauge fixing condition (2) Lagrangian (1) has the form: 

C = ^A m d ll A m - ^A AA - Ao]o + A m j m - ^{i-fd^ - m)^ 
Substitute into this Lagrangian the solution of the classical equation of motion (3): 

S = J d 4 x - B.Bi) + Aiji - ^7^ - m)il> + ^jo^Jo) (4) 
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3. POSITRONIUM EFFECTIVE ACTION 



The quantization of fermion fields in action (4) yields generating functional: 



Z = Jti \ D^D^e J , (5) 



Ai here and below we consider only as external fields. 
In this generating functional consider the term in (4) 

i J d 4 x d 4 y j^(x)D^(x - y)j„(y) = Jd 4 x ^ j ^j 



—Jd^d^d^d^ ^ ai (x 1 )^(x 2 ) 

7 a 2 5 W-X2) 



47r|xi-y 2 



5 3 (yi-y 2 )7 g 



/c " 1 /3 1 o 2 ' 32 (zi,yi;z2,y2) 



^ 2 (xly 2 )^(x° 1 ,y 1 ) = 



where in the first line D^ u (x — y) — —Sn5n ; :— , and in the last line we use next 

y 1 °° 47r| Xl -y 2 | 

notation: sum over pair indexes a±,bi is summing over ai, fli and integrating over dfaidfyi. 

So, there is forth power of ip in generating functional. Habbard-Stratonovich transforma- 
tion allows to reduce power of ip to second. 

Let's consider the combination ipip as a real bilocal field: 

X yx°,x,y) = V Q (Ax)^(x°,y) 

Transposition operation we define as: 

t fit o \ i o \ 

X a {x ,x,y) = -X a {x ,y,x) 



Than introduce new bilocal field M a /3 (x°, x, y) and make 
1 . 1 



i Id x - jo— jo 



2 J A J = e 2 



V det /C 1 is constant and has been included in ^/K. 
After that generating functional (5) takes form 



Z = J/ 2 / Dt/'Dt/'DjM exp i 



<l L .r [ ^{A t Ai ~ BiBi) + Aiji - cd-i'i),, - m)c j + 



+\M T ai bl IC- ia \ ia2 b2 M a \ 2 +^ b M\ + J d*x (ftf + jrj) 



Introduce notation: 



G m \%(x, y) = (-<y + mb % + - V) + M a p(x°, x, y)5(x° - y°) 

and define inverse operator as: 

/ dSjG-\ a p (x, y)G m / 1 (y, z) = 5%5\x - z) 
that allows us to write generating functional in more compact form: 



Z = J/ 2 DipDipDM exp i 



M a \+ 



After integration over fermions finally we have: 



Z = J/ z / DM exp i 



-M T bl r- iai 62 M a2 - 



J d 4 xd\jr} a (x)G mA a '- p {x, y)r] p {y) - itr In G m \ 



So, effective action for positronium field take form: 
S P 



= J d*x Q(M - B^)) + i-M T ai 6l /C- iai 6ia2 62 ^- b2 - itr In 



(6) 
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4. QUANTIZATION OF BILOCAL FIELDS 

For our purpose it would be enough to use the semi-classic approach. According to this 
approach the first variation of the action (6) is the Schwinger-Dyson equation : 

that gives a fermion spectrum. And the second variation in the point of minimum is the 
Bethe-Salpeter equation: 



5 2 S 



p 

r = o 



SM 2 

which allows us to find a bound state spectrum 



4.1. Schwinger-Dyson equation 

Denote solution of (7) x, y) — m5 a /3 (5(x — y). Than introduce notation: 

G?%(x,y) = -i-y^fa + S%(x ,x,y)5(x° - y°), 
where the inverse operator is defined as 

J d\Gz ia p(x,y)Gx%(y,z) = 8%8\x - z). 
The equation (7) now takes the form: 

Z a \ 1 =mS\+iJC a \ ia2 b >Gx a \ 2 , (8) 

in the last term the sum is over pair indexes 02,62 is summing over 02, /?2 and integrating 
over d 4 x 2 d 4 y28(x2 — y%)- Below we also assume a delta-function S(x® — y°) m contraction of 
/C and G. 

After Fourier transform let us takes an ansatz: 

Z(p, q ) = 5\p-q) (q + £(q)e- 2 ^ (q) ) . 
Then we can inverse G^ 1 : 

f A+(q) A-(q) \ (q) 

S \ qo -E(q)+te ^ qo + E{c i )-ie) 1 ' {) 



where: 



A ± (q) = e A^)l±l! e -4^). (10) 
2 
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Well known that in case of QED Schwinger-Dyson equation gives small correction to 
electron mass, so we take the solution of Schwinger-Dyson equation in form S a 6 mS a b , 
where m is the electron mass. 

4.2. Bethe-Salpeter equation 

Taking the second variation of the action (6) in the point of minimum that we have found 
solving the Schwinger-Dyson equation (8): 

S 2 S f 



5M 2 

we obtain the Bethe-Salpeter equation: 



r = 

A=0,M=T.-m 



r ai , = iJC a \ n b2 G^ a \ T a \ G^ b \ . 

oi oi a2 ^ a 3 03 ^ 02 

After Fourier transform and inserting Gs in (9), let us take integral over q : 

r(p) = J d 3 q/C(p-q)^(q) (11) 

where the wave-function: 

» (q ) = V P + (l)r(q)A-(q) + 7C-(q)r(q)A + (q) \ (12) 

VH; ' V Ep-mp + ie E P + m P - ie J ' y ' 

where analogously to (10) one has: 

A ± (q) = e -^ ( ^e4^. (13) 

Acting (13) and (10) on (12), and inserting (11) in (12) we have equation for \E': 

(£ pT m P )A ± (p)M/(p)A T ( P ) = A±(p) (J d 3 q/C(p-q)^(q)^) A T (p) 

If |p| <C m than this equation turns into the Schrodinger one [12]: 

^p 2 + (2m - m P ^j *(p) = J d 3 q/C(p-q)M>(q) (14) 

As we interested in parapositronium, we can take the next ansatz: 

H>(q) = iV 7 5 7V(p) (15) 
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where iV is unknown normalized factor and will be defined below, ^(p) is the nonrelativistic 
wave-function - the solution of Schrodinger equation (14): 



m«\ 2 _ r rna 

e 2 

7T V 2 / 



^(P) 



1 



(ma) 2 



(16) 



'((?)' + pf 

Because of smallness a the condition |p| <C m is satisfied. Therefore the motion of an 
electron and a positron in the positronium is nonrelativistic and we can take Mp ~ 2m. 



5. TRIANGLE ANOMALY 



Small variation action Sp (6) over Ai and double variation over A give us 



M'A m A k 



5 5 



SM 5A m 5A k 



S, 



-le 



A=0,M=^-m 



that can be interpreted as the triangle anomaly (fig.l) with an amplitude 



M, 



dq tr 



G s (g + |) r(g) G E (g - |) i x (* x ) G S (^A + g) A 2 (fc 2 ) 




Below in our calculations we consider one virtual photon with momentum k% and other 
real photon with momentum k,2- For simplicity the reference frame is chosen as a positronium 
rest frame V = (Mp, 0, 0, 0). Let us denote k 2 as k. 

Rather crude estimate can be done in the following way. Using (9), let us integrate over 
go and take into account only a pole that forms the wave-function (12), using (15) 



/ 



M A = e 2 N / ^q^(q)tr 



k 2 - h 



7 5 7°Ai(fci 



+ q + m 



\ 



V 



2 ( |k| a/ m 2 + q 2 — kq 



A 2 (A;2 
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Let us make a change 7 for ^ and make the trace over the gamma matrixes product 



M A = -Ae 2 N J d 3 q^(q)e 0Xpu e»e 



Z [ 

2 2 (jk|^/m 2 + q 2 -kq) M ^ ' 



The term with g A gives zero contribution. We can change Eo\^ u —— for sum e p \ pv —- , 

Mp Mp 

furthermore we can suppose now that in arbitrary frame of reference (not in positronium 
rest frame of reference as we do before) the last sum should appear. Using (16) one gets: 

M A = Ae 2 N e, uXp k»e»k& / d\ \ - 2 - -J— . 

Notice, that as a is small, so only values of q near zero are important. Therefore, we can 
write in spherical variables, where £ = — cos 9, 

00 1 

f f q 2 1 

2tt / dq / dx 



J 7 ^nu*y + q2 y (ky/m 2 + q 2 - kqC) 



2tt J dx —. ^ = ^ = r- f 



> x (k^Jm 2 + q 2 — A;gC j 

and 



oc 



rfg g = 2?r- 



((^f) 2 + 4 2 )' 



km 2am' 



_ 27Va§ Fe ,,A e P 



To get agreement with the decay P — > 77 given by anomaly, we should take the normalization 
factor as 

Thus finally we have for the triangle diagram 

M A = C P e^ Xp k?e» 1 k$e p 2 , (17) 

where 

Cp ~ Am 2 - k\ ~ (k 2 V) ■ (18) 
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6. HAMILTONIAN ANOMALIES 

In quantization of the Hamiltonian systems we need canonical coordinates and momenta. 
If the triangle anomaly exists, it changes the momentum. Consider as an example: 

which is the action (6) only with the first term and the triangle anomaly (17) with V = 
(Mp, 0, 0, 0) and Cp = C - consider as a constant, <fi is the wave function of positronium. 
The canonicals momentum: 

so, one can write: 

and consider the first term as a "renormalized" electromagnetic field and the second term 
as a Hamiltonian anomaly shown in fig. 2. 



A'i 




Figure 2. Hamiltonian anomaly. 

Back to the action (6) and triangle anomaly (17,18), the amplitude of the Hamiltonian 
anomaly (fig. 2) is 
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Evaluating sum over v 

lCp(p 1 ,k 1 )C P (p2,k 2 ) , x 
M x = — 2 {Ve 1 p e 2p + T» e^e 2v ) , (19) 

where 

V = -(pip2)(fcifc2) + {pik 2 )(kip 2 ), ^ 
= (piP^WkS - {k x p 2 )k 2 » Vl v - {p x k 2 )p 2 »k x v + {k x k 2 )p 2 »p x v . 
The square of the amplitude equals to 

2 n 2 a 10 ( 

l M x| = ~77, 7T7-, T2 2 (PlP2)(fclfc2)bl^l)(P2^2) +2(piA:2)(^lP2)(Pl^l)(P2^2)- 

4(fcipi) (A; 2 p 2 ) V 

- 2M P 2 (p 2 k 1 ){p 2 k 2 ){k 1 k 2 ) - 2Mp 2 (k 1 k 2 ){p 1 k 1 ){ Pl k 2 ) + M P A {k l k 2 ) 

After an integration over the positronium momenta p\ and p 2 , finally the cross-section is 



(2s + 7M P 2 )Wl- -6M P 2 ln 



where s — (k\ + & 2 ) 2 and Mp - positronium mass. 

Consider the behavior of a at small and large s in the center of mass frame of reference. 
Denote the energy of one photon as E 1 , then s = 4E" 7 2 . For E 1 —¥ M P we have: 



7ra 10 i fM P ^ 2 



128M P 2 V 1 V E~ 1 



j 1 



and for E~ — > oo we have: 

cr — >■ 

48Mp 2 

In conclusion we should notice that in the process 77 — > PP, where P is the parapositro- 
nium besides (19,20) there are leading contributions from the box diagram. 



7. CONCLUSION 

In this paper we obtain the bound state generation functional by Poincare-invariant irre- 
ducible representations. Specific feature of this method is an absence of explicit relativistic 
covariance. An attempt to conserve the covariance based on Markov- Yukawa constraint was 
considered in [13]. We will prolong this research in future papers. Some comments about 
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the relativistic covariant excluding the time component discussed in section 2 can be found 
in the review [6]. 

The quantization of the positronium action (6) is made in the semi-classical approach. 

As an example of this method, we calculated the triangle anomaly (fig.l). To simplify the 
calculations the assumption about poles was made. Besides in resulting the answer (17,18) 
was made an assumption about relativistic structure of the answer. 

The anomaly triangle diagrams change canonical momenta in Hamiltonian systems. That 
may leads to additional contribution (19,20) to the process 77 — > PP. We can calculate this 
contribution using the previous results with the triangle diagram. 
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